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Linear  functional-differential  equations  >;ith  constant  coefficients. 


Jack  K.  Bale 


I,  Introduction.  Until  recently,  most  of  the  results  concerning 
differential-difference  equations  have  been  obtained  by  treating  the 
dependent  variable  as  a  point  in  Euclidean  space  and  employing  arguments 
which  are  standard  in  the  theory  of  ordinary  differential  equations* 

To  the  author fs  knowledge,  Krascvskii  [  9  ]  was  the  first  to  exploit 
the  idea  that  the  proper  setting  for  these  problems  is  in  a  function 
space.  In  doing  so,  the  argunents  used  for  ordinary  differential  equa¬ 
tions  become  more  natural  for  differential-difference  equations.  The 
present  paper  is  an  attempt  to  obtain  some  analogies  between  linear 
differential-difference  equations  with  constant  coefficients  and 
ordinary  linear  differential  equations  with  constant  coefficients. 

Mon  specifically,  we  discuss  in  detail  the  eige  ns  paces  of  the 
linear  equation  and  then  make  use  of  the  adjoint  equation  to  introduce 
new  coordinates  In  the  function  space  which  exhibit  In  a  natural  manner 
the  behavior  of  the  solutions on  an  eige ns pace  and  the  behavior  of  the 
solutions  on  a  complementary  space.  In  this  manner,  it  is  shown  in 
section  IV  how  many  of  the  usual  perturbation  theorems  in  ordinary 
differential  equations  can  be  easily  extended  to  differential-difference 
equations.  The  basic  idea  for  discussing  the  problems  in  this  manner  is 
contained  in  the  papers  of  Shimanov  (see  the  bibliography)  and  the  present 
paper  originated  from  an  attempt  to  understand  the  geometric  significance 
of  Shimanov's  results.  This  approach  should  lead  to  a  better  understanding 
of  much  of  the  geometric  theory  of  differential-difference  equations.  The 
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author  is  indebted  to  John  Steulpnagel  and  Arnold  Stokes  for  many 
fruitful  discussions. 

The  following  notation  will  be  used  throughout  this  paper.  Rn 

is  the  linear  space  of  n-vectors  and  for  x  €  Rn,  \y:\  is  any  vector 

norm.  For  any  given  numbers  a,  0,  a  £  p,  C(  [a,  pi,  Rn)  will  denote 

the  space  of  continuous  functions  napping  the  interval  [a,  0]  into  Rn 

and  for  9  €  C(  [a,  p],  Rn),  \\tp\\  *  sup  |q>(0)|.  For  any  r  2  0, 

a  5  9  £  0 

any  continuous  function  x(u)  defined  on  -r  £  u  £  A,  A  >  0,  and  any 
fixed  t,  0  £  t  £  A,  we  shall  let  the  symbol  x^  denote  the  function 
xt(0)  *  x(t  +  9),  -r  £  0  £  0;  that  is,  x^  €  C(  [a,  0],  Rn)  and  is 
that  "segment”  of  the  function  x(u)  defined  by  letting  u  range  in 
the  interval  t  -  r  £  u  £  t. 

Let  X(q>,  t)  €  Rn  be  a  function  defined  for  all  9  €  C([-r,Q],  Rn), 
II9II  £  H,  H  >  0,  t  €  [0,»).  Let  x(t)  denote  the  right  hand  derivative 
of  a  function  x(u)  at  u  =  t,  and  consider  the  functional-differential 
equation 


(1.1)  x(t)  =*  X(xt,  t). 

Definition  1.1.  let  tQ  be  any  given  number  k  0  and  let 
9  €  C( [-r,  0],  Rn),  ||9||  £  H,  be  any  given  function.  A  function  x.(t  ,  9) 
is  said  to  be  a  solution  of  (l.l)  with  initial  function  9  at  t  if  there 
is  a  number  A  >  0  such  that 


i)  for  each  t,  tQ  £  t  £  tQ  ♦  A,  \(^0>  9)  is  defined, 
belongs  to  C([-r,  0],  Rn)  and  ||xt(to,  q>)||  s  H; 
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ii)  (*0>  <p)  -  <p; 

o 

lii)  x(tQ,  9)  satisfies  (1.1)  for  tQ  £  t  £  tQ  ♦  A. 

If  t  is  equal  to  zero,  we  shall  abbreviate  x(tQ,  9)  by 
x(9).  If  X(9>  t)  is  continuous  in  9,  t  and  Lipschitzlan  in  9,  it 
is  easy  to  prove  that  (1.1)  always  has  a  solution  and  for  each  9  there 
is  only  one  solution.  Furthermore,  it  is  also  easy  to  prove  that 
x(tQ,  9)  Spends  continuously  on  9. 

By  a  linear  functional  ^differential  equation  with  constant  coeffi¬ 
cients,  we  mean  a  system  (1.1)  where  X(q>,  t)  =  f(q>)  is  homogeneous  and 
additive  in  9.  It  is  well  known  [15,  p.  110]  that  f(qp)  continuous  cm 
C([-r,  0],  Rn)  implies  there  is  a  matrix  tj(0)  whose  elements  are  of 
bounded  variation  such  that 

0 

f(<p)  =  /  [dn(e)  to(e), 

-r 

for  all  9  €  C(  [-r,  0],  Rn),  where  the  integral  is  in  the  sense  of  Stieltjes. 
This  observation  makes  it  obvious  that  the  concept  of  linear  functional- 
differential  equation  with  constant  coefficients  includes  all  linear 
differential-difference  equations  with  constant  coefficients  of  the  form 


±(t) 


P 

*  L  Ak  x(t  -  T^), 
k*l 


2  0. 


-  k  - 


II*  Basic  properties  of  linear  systems  with  constant  coefficients,  A 
linear  functional-differential  equation  with  constant  coefficients  is 
any  equation  of  the  form 

(2.1)  i(t)  =  f(xt) 

where  f  is  a  continuous  linear  function  mapping  C(  [-r,  0  ],  Rn)  into 
Rn.  For  any  such  function  f(q>),  it  is  well  known  (see  [15,  p.  110]) 
that  there  exists  an  n  x  n  matrix  i\(0),  -  r  £  0  £  0r  whose  elements 
have  bounded  variation  such  that 

0 

(2.2)  *(<p)  -  /  [dn(e)Me). 

-r 

If  9  Is  any  given  function  in  C( [~r,  o],  Rn)  and  x(q>)  is  the 
solution  of  (2.1)  with  the  initial  function  9  at  zero,  we  define  the 
operator  Cf(t)  mapping  C(  [-r,  0],  Rn)  into  C([-r,  0],  Rn)  by  the 
relation 

(2.5)  xt(<r)  =  O^(t)q), 

where,  for  each  fixed  t  £  0,  ^(9)  is  the  function  in  C([-r,  0],  Rn) 

determined  by  the  relation  x^(y)(e)  -  x(9)(t  +  0),  -  r  £  0  £  0. 

Lemma  II.  1.  The  operator  ^T(t),  t  i  0,  defined  on  C([-r,  0],  Rn) 
by  (2.3)  satisfies  the  following  properties. 
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i)  7  (t)  is  a  bounded,  linear  operator  for  each  t  2  Oj 

ii)  Sf  (t)  is  strongly  continuous  on  [0,  »);  that  is 

,7(0)  =  I  and 

liJ\  _» J  CTCO*  -  Jr('t)q> II  =  o, 

for  all  t  i  0,  <p  e  C([-r,  0],  Rn); 

iii)  The  family  of  transformations  (  J"(t),  t  i  0)  Is  a  semigroup, 
that  is, 

»7( t  +  T)  =J(t)7(x)i  a11  t  i  0,  T  2  0; 

iv)  J  (t)  is  completely  continuous  (compact)  for  t  i  r;  that 

Is  ,  J(  t),  t  £  r  maps  closed  bounded  sets  into  compact  sets. 

Proof:  i)  It  is  obvious  that  ^(t)  is  linear.  Since  f(q>)  is  con¬ 
tinuous  and  linear,  it  follows  that  there  is  a  constant  L  such  that 
jf (9)  |  £  L||q>||  for  all  From  the  definition  of  *7(t),  we  have,  for 
any  fixed  t 

%7  (*)$(&)  =  <p(t  +  9),  t  +  $  £  0, 

(2.4)  t+e 

J  (t)q>(e)  =  $>(0)  +  /  f(  J  (t)<p)dT,  t  +  9  >  0,  -r  M  <  0. 

0 

Since  |f (qp)  |  i  L\\yW,  it  follows  that 

II  7(t)<p||  5  «Ulk>ll,  tiO,  9  e  C([-r,  0],  Rn), 
and  7(  t)  is  bounded. 
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ii)  From  i),  it  follows  that  ,/(t)  is  continuous  for  all 

t  20  and,  from  the  definition  of  3  (t),  3( 0)  «  I.  This  proves  ii). 

iii)  This  is  immediate  from  the  definition. 

iv)  To  prove  iv),  we  observe  that  if  S  »  {9  €  C([-r,  0],  Rn)  |  \\q\\  S  R) 

then 

'7  (t)sC  B1  -  It  €  C([-r,  0],  Rn)  Ii  €  C([-r,0],Rn),  M  S  eLtR, 

lllll  s  I*LtR} 

for  t  2  r.  Since  is  compact  and  J (t)  is  continuous,  the  result 
follows.  This  completes  the  proof  of  Lemma  II.  1. 

For  any  semigroup  of  transformations  3( t),  t  *  0,  of  a  Banach 
space  into  itself,  the  infinitesimal  generator  ^  of  t)  is  de¬ 
fined  by  the  relation 


=  +  ^(t)<p  -  <p] 

t  ->  0 


for  every  value  of  9  for  which  this  limit  exists.  The  limit  of  course 
signifies  convergence  in  the  norm  of  10  . 

For  any  operator  3  of  a  Banach  space  >0  Into  itself,  the  resolvent 
set  p(J)  of  3  is  the  set  of  values  \  in  the  complex  plane  for  which 
the  operator  XI  -  3  has  an  inverse  which  is  defined  for  all  9  in  . 
The  complement  of  p(  {/)  in  the  complex  plane  is  called  the  spectrum  of 
and  is  denoted  by  v(3)*  The  spectrum  v{3)  of  an  operator  consists 
of  three  different  types  of  points,  namely  the  residual  spectrum  K<j(3)> 
the  continuous  spectrum  Ca(3)f  and  the  point  spectrum  Pa (3)*  The 
residual  spectrum  consists  of  those  values  of  X  in  a(^)  for  which 


XI  -  3  exists  but  the  domain  A?(xi  -  %f  )-1  of  (xi  -  3  )-1  is  not 
dense  in  ^  .  The  continuous  spectrum  consists  of  those  X  in  a(  ) 
for  vhich  £(\I  -3 )-1  is  dense  in  35  and  xhe  point  spectrum  consists 
of  those  values  of  X  for  XI  -  3  does  not  bare  an  Inverse.  The  points 
X  In  P a(3 )  are  sometimes  called  the  eigenvalues  of  £f  and  any  non¬ 
zero  9  such  that  (XI  -£f  )q>  =  0  is  called  an  eigenvectors 

One  of  our  first  objectives  is  to  try  to  determine  the  imture  of 
a(J(t))  and  o(^)  for  the  family  of  operators  vhich  arise  in  our 
particular  problem  and  to  analyze  in  what  sense  the  operator  Cf  (t)  is 
approximated  by  the  operator  e^*  provided  this  latter  object  makes 
sense.  For  the  simple  case  in  vhich  system  (l)  is  an  ordinary  differential 
equation;  that  is,  f(<j?)  =  A9(0)  for  some  constant  matrix  A,  the 
operator  (t)  is  eA^  and  the  infinitesimal  generator  of  3 (t) 

is  equal  to  the  matrix  A.  The  following  results  show  that  analogous 
results  are  valid  for  the  more  general  system  (2.1). 

The  following  lemma  is  a  restatement  of  Theorem  10.3.1  and  10.3*3  of 
HiUe  and  Phillips  [  8  ]  for  out  particular  case. 


lemma  II. 2.  If  J(t),  t  2  0,  is  a  strongly  continuous  semigroup  of 
operators  mapping  C([-r,  0],  Rn)  into  C(  [-r,  0],  Rn),  then  the  domain 
of  the  infinitesimal  generator  ^  of  J^(t)  is  dense  in 
C([-r,  0],  Rn)  and  the  range  of  $  is  in  C([-r,  0],  Rn).  For 

all  9  in 


(2.5) 


jjr  3 (t)9  *  3 (t)^9  “ 
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We  now  derive  a  specific  formula  for  the  infinitesimal  generator  C 
in  terms  of  the  system  (2.1).  Since  x(9)  satisfies  (2.1),  it  follows 
from  the  definition  of  £f ( t )  that  Cf  { t)9  satisfies  relation  (2.4). 

Consequently,  for  any  0,  -r  £  e  <  0, 

U*  +  9(0)]  =  3ggl. 

t  0 


If  0  =  0,  then 


11m  +  r[(<7?t)<p)(0) 

t  -»0 


9(0)1-  lim.k/  ffrjdr] 
t  -»0  0 


-  f(xQ)  =  f(<p) 


=  /  tdr,(e)]9(e). 


-r 


But,  to  say  that 


=  lim  .  r[<^(t)<p  -  9] 
t  ->0 


implies  convergence  in  the 
convergence  and  thus  9 
^<p(0)  =  f(9)  =  d9(O)/d0. 
lng  lemma. 


norm  in  C([-r,  0],  n),  which  implies  uniform 

must  he  a  continuous  function  which  implies 
Summarizing  these  remarks,  we  have  the  follcrw- 


Lemaa  II. 3.  If  v/(t),  t  20  is  the  family  of  transformations  on 
C([«r,  0],  Rn)  defined  by  (2.3),  then  the  infinitesimal  generator  of 
ijoo)  is  given  by 


9(0)  = 


^  -hseso 

;°r[dTl(0)]9(0)  =  f(9). 
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tie  now  take  some  results  from  Hills  and  Phillips  and  Reiez-Hagy 
and  apply  them  to  our  family  of  operators  ^/(t)  satisfying  properties 
i)  -  It)  of  Lemma  II. 1  to  obtain  the  following  results. 

Lemma  II. 4.  For  t  2  r,  o{\J( t))  is  a  countable  set  and  is  a  compact 
set  of  the  complex  plane.  The  only  possible  point  of  accumulation  of 
a(jf(t)),  is  (0}  and  if  fi  /  0,  4  in  o(J( t)),  then  n  is  in 
Pa(^(t)).  (Hi lie  and  Phillips  [8,  pp.  180-182],  Theorems  5*7*1,  5-7-3]) • 

Given  any  operator  0  ,  we  denote  by  {£])  the  null  space  of 
Cf  }  that  is,  the  set  of  all  9  such  that  2^ 9  =  0. 

Lemma  II. 5.  For  t  k  r,  if  \i  =  *i(t)  is  in  Po(^ (t)),  \i  /  0,  then  for 
each  positive  integer  k,  (|j.I  -  v7(t)k  is  of  finite  dimension  for 
every  k,  and  there  exists  a  least  integer  nQ  such  that 

utte  -cr(t))k=  for  all  k,  |  *  nQ. 

If  W(J(t))  is  equal  to  ^(ni  -v7( t))n°,  then 

(Hille  and  Phillips  [8,  p.  182,  Theorem  5*7*3]) 

Lemma  II. 6.  For  t  £  r,  Pc(  \f  (t))  =  exp[tPa(<^)]  plus  possibly  {0). 
More  specifically,  if  4  =  y(t)  is  in  Pa(^/(t) )  for  some  fixed  t  and 
H  /  0,  then  there  is  a  point  X  in  ?a(CL)  such  that  Further¬ 

more,  if  (Xn)  consists  of  all  distinct  points  in  Po^O  such  that 
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e  n  =*  y,  then  fit.  (4I  -  is  the  linear  extension  of  the  linear 

independent  manifolds  (  ^(X^I  -  £?)*).  (Hille  and  Phillips  [8^  pr  467, 
Theorem  16.7.2  and  pp.  521-524,  Theorem  10.6. 5]). 

In  the  following,  ^(B)  for  any  operator  B  on  a  Banach  space 
and  X  in  Po(B)  will  denote  the  maximal  subspace  of  ^  annihilated  by 
powers  of  B  -  XI. 

Lemma  H.7.  If  ^  is  the  infinitesimal  generator  of  the  family  of  operators 
defined  by  (2.5)  and  X  is  in  P a(^ ),  then  the  set  in 

C(  [-r,  0],  Rn)  is  finite  dimensional.  Furthermore,  there  is  a  real  number 
0  such  that  Re(X)  £  0  for  all  X  in  Pa(^),  and  there  are  a  finite 
number  of  X  in  such  that  r  *  Re  X  for  any  given  real  number 

r. 

This  lemma  is  an  Immediate  consequence  of  Lammas  II. 4  -  II. 6. 

Lemma  II. 6  above  gives  a  very  distinctive  relationship  between  the 
point  spectrum  of  Cf  (t)  and  the  point  spectrum  of  ^  .  In  fact,  except 
for  the  point  \t  =  0,  Pc(v7(t))  is  completely  determined  by  Pa(^).  We 
now  derive  an  explicit  expression  for  Po(^)*  If  X  is  in  Pa(^),  then 
there  must  exist  a  nonzero  9  in  C( [-r,  0],  Rn)  such  that 


X  ?  =  xp* 


This  last  relation  is  satisfied  if  and  only  if 

=  Xflp(0),  -r  i  8  s  0 
1°  Mo)  Me), 
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which  in  turn  lns>lles  9(0)  =  e*  b,  -r  S  0  i  0,  and  the  rector  b 
satisfies 


(2.6) 


A(x)b  =  0 


0 


A(X)  =  (XI  -  /  (dri(0)]eX0). 

-r 


As  a  result  of  this  fact  and  Lemma  II. 7,  we  have 


Lemma  n. 8.  P a(fj)  ~  {x|  det  A(x)  »  0}.  The  roots  of  the  characteristic 
equation  of  (2.1),  det  A(x)  =  0,  have  real  parts  bounded  above  and  there 
are  only  a  finite  number  with  real  parts  greater  than  a  given  constant. 

The  characteristic  equation  of  (2.1)  can  be  obtained  in  a  very 
straightforward  manner  as  in  ordinary  differential  equations  by  determining 

necessary  and  sufficient  conditions  that  the  equation  (2.1)  has  a  solution 

/  >.  Xt 

of  the  form  x(t)  s  e  b  for  some  X,  b.  Many  procedures  have  been  given 
to  analyze  the  nature  of  the  roots  of  the  characteristic  equation  (see, 
for  example,  Ianger  [11],  Pontrjagin  [14],  Pinney  [13]).  The  property 
of  the  roots  mentioned  in  Lemma  II. 3  is  also  in  these  papers. 

Before  proceeding  to  an  analysis  of  the  specific  structure  of  the 
solutions  of  (2.1),  we  wish  to  derive  a  result  concerning  the  maximum 
rate  of  growth  of  the  solutions  of  (2.1).  For  later  reference,  we  state 
the  result  independently  of  solutions  of  (2.1).  If  is  a  bounded  linear 
transformation  of  a  Banach  space  into  itself,  the  spectral  radius  p  y  of 
is  the  smallest  closed  disk  with  center  at  the  origin  in  the  complex 
plane  which  contains  a(J ).  The  following  lemma  is  taken  from  Reisz-Nagy 
[15,  P.  425]. 
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Lemma  11,9*  If 


3  u 


a  bounded  linear  transformation  of  a  Banach  space 


into  itself,  then  the  spectral  radius  p j  of  is  given  by 


=  lim  ||Jnj|1/n. 

n  -+  » 


Theorem  II.  1.  If  Cfi  t), 
operators  of  a  Banach  space 


t  2  0,  is  a  strongly  continuous  semigroup  of 
nS  into  itself,  if  for  some  r  >  0,  the 


spectral  radius  p  *  ^(r)  is  ai*d  t  0  and  0r  *  log  p,  then,  for 

any  y  >  0,  there  is  a  constant  K(y)  fe  1  such  that 


||^(t)9||  ^  K(r)e^^r^|j9||  for  all  t  JO,  qp  In  10* 


Proof:  This  proof  is  essentially  the  same  one  as  contained  in  Stokes 
[19].  Since  p  is  finite,  the  number  p  is  veil  defined.  Prom  Lemma  II. 9, 

=  11m  i|Jn(r)||1/n,  and  thus,  for  any  y  >  0, 

n  -» » 

e-^  =  lim  e-(p+r>r||JIi(r)||1/n. 
n  -4® 


Therefore,  there  exists  a  number  N  such  that 


e.(Mr)nr||(7»(r)| 


+  €  ) 
n' 


n 


where  e-TT  +  €n  ^  K  <  1  for  all  nU.  Consequently,  e^^^nr\\^fn(r)\\  -4  0 
as  n  -►  ®.  Since  is  continuous  for  all  t  fc  0,  there  is  a  con- 

stant  B  ouch  that  ||J(t)||  SB  for  0  S  t  s  r.  Define  K(r)  for  any 


y  >  0  to  he 
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K(r)  -  ■*{  nax  Be-(P+r)nr||  Jn(r)ll). 

0  £  t  *  r  n  fc  0 

If  0  £  t  £  r,  then,  for  any  9  in  ^ , 

li^(t)»||  *  ll^(t)||-|MI  *  B||<p||  S  K(r)e(p+r)tll9l|. 

If  t  £  r,  then  there  is  an  integer  n  such  that  nr  S  t  <  (n+lr), 
and,  for  all  9  in  , 

II  C7(t)q>i|  =  ||#(t  -  nr)(f  (nr)<p||  5  B||(*/n(r)||»  |jcp||  s 

[Be-(^)(t^)e-((W)«ir||yn(r)y]e(p*r)tw 

S  K(r)e(&+T)t||9||. 

This  completes  the  proof  of  the  theorem. 

noyollary  TT.l.  if  A(X)  is  defined  as  in  (2.7)  and  all  the  roots  of 
the  characteristic  equation  det  A(x)  =  0  satisfy  Re  X  £  0,  then  for 
any  y  S  0,  there  exists  a  constant  K(y)  &  i  such  that  if  x(9)  is  the 
solution  of  (2.1)  with  initial  function  9  in  C(  [-r,  0],  RU)  at  zero, 
then 


ll*t(<P)ll  *  K(r)e(^)t|h>l|,  t  *  0. 

In  particular,  if  0  <  0,  then  all  solutions  of  (2.1)  approach  zero 
exponentially  as  t  -s>«. 

Proof:  If  47 (t)  is  defined  as  in  (2.3),  then  p  «  p finite  and  *  exp  0r 
from  Lemmas  II. 6  and  II. 8  if  p  /  0.  This  case  follows  from  Theorem  II.  1. 

If  p  =  0,  the  corollary  is  obviously  true. 
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Corollary  II.  1  Is  veil  known  in  the  literature  (see  for  ex&ngjle, 
Krasovski  I  [9],  Bellman  and  Cooke  [  1],  Stokes  [19]). 

In  the  following,  we  shall  always  assume  that  \^(t)  *s  family 
of  operators  associated  with  system  (2.1)  and  ^  is  the  infinitesimal 
generator  of  t ) . 

From  Lemma  II. 7,  if  X  is  in  Po£^),  then  <5^ (XI  -^0*  is  of 
finite  dimension  for  every  integer  k  and  there  is  an  integer  nQ  such 
that  (XI  -^)k  =  J^(XI  -fl)  °  for  all  k  £  nQ.  Since 

-  M)  °  has  finite  dimension,  say  d,  there  exists  a 
basis  q»x,  9d  of  Furthermore,  since  Q 

there  exists  a  d  x  d  constant  matrix  B  (whose  only  eigenvalue  is  X) 
such  that 

=  4>B,  $  =  row^,  fd). 

Now,  let  us  investigate  the  nature  of  the  solutions  of  (2.1)  given  by 
(t)fl>.  From  Lemma  II. 2, 

^  =  J(t)4»  =  J(t)*B 

and,  thus, 

=  le®1. 

Consequently,  we  have  the  following  result. 

Theorem  II. 2.  Let  -7(t),  t  £  0,  be  the  strongly  continuous  semigroup 
of  operators  defined  by  (2. 3)  and  let  {7  be  the  infinitesimal  generator 
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of  t).  Suppose  X  is  in  P o(^)  and  let  4>  =  rowfo^  ^d) 

tie  a  basis  for  w^d),  the  maximal  subspace  of  C([-r,0],  Rn)  anni¬ 
hilated  by  powers  of  XI  -  If  a  is  an  arbitrary  constant  colunn 
vector  of  dimension  d,  then  J(t  )$a  =  ♦e^a  for  some  constant  matrix 
B  with  all  of  its  eigenvalues  equal  to  X.  If  n  *  *i(t)  ^  0  is  in 
Pa(J(t)  for  some  t,  and  (Xp  . ..,  Xp}  is  the  set  of  distinct  ele¬ 
ments  in  Pa(^f)  such  that  e^  =  n,  and  $  =  (*p  . ..,  *^),  where 
is  a  basis  for  (^)>  <iiioenflion  d^  then,  for  any 

<P  €  ty(Cf{t)),  there  exists  a  vector  a  of  dimension  d  «  +  •  •  •  +  d^ 

such  that  q>  =  #a  and 


(2.7) 


=  «  eKa 


where  B  is  a  dxd  matrix  given  by  B  =  diag  (Bp  . ..,  B^)  where 

each  B.  is  defined  by 

J  j  J  J 

The  first  part  of  this  theorem  has  been  proved  above  and  the  second 
part  follows  by  using  the  same  argument  together  with  Lemma  II, 6. 

Theorem  II. 2  shows  that  on  the  subspace  <$^(7(t)),  the  functional- 
differential  equation  (2.1)  has  the  same  structure  as  an  ordinary  linear 
differential  with  constant  coefficients.  Of  course,  the  dimension  of 
the  matrix  B  in  (2.7)  may  have  nothing  to  do  with  the  dimension  of 
syBtem  (2.1)  since  the  dimension  of  B  is  determined  by  the  multiplicity 
of  the  eigenvalue  ^(t)  of  Notice  also  that  the  multiplicity  of 

this  eigenvalue  may  change  with  t.  In  fact,  this  is  easily  illustrated 
by  an  ordinary  equation 

2rri  0 

i«Ax,  A-(0_arl). 
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Then  v7(t)  ■  and  Pa(\7 (t))  «  (e^'S  e-£)rtt^  For  t  »  n  an 
Integer,  Po(\7(n)  *  (1}  and,  otherwise,  consists  of  tvo  distinct  points. 
On  the  other  hand,  the  multiplicity  of  the  points  in  Po(  (?()  do  not 
change  and  one  can  always  define  the  set  to  generate  solutions 

of  (2.1). 
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HI.  The  adjoint  equation  and  a  change  of  coordinates.  In  this  section, 
we  wish  to  show  how  one  can  introduce  a  change  of  coordinates  in  the 
space  C(  [-r,  0],  Rn)  in  such  a  way  as  to  exhibit  in  a  natural  manner 
any  particular  eigenspace  Wjd)  associated  with  system  (2.1).  More 
specifically,  we  show  that  it  is  possible  to  transform  system  (2.1)  to 
an  equivalent  system  which  consists  of  a  set  of  ordinary  differential 
equations  (whose  solutions  describe  the  behavior  of  //(t)  on  the  eigen¬ 
space  KTfy(Sjl))  together  with  an  operator  equation  (whose  solutions 
describe  the  behavior  of  SJ (t)  on  a  space  complementary  to  WaH)). 

To  do  this,  we  make  use  of  the  equation  "adjoint*1  to  (2.1).  This  concept 
has  been  used  in  functional -differential  equations  by  many  authors  (see, 
for  example,  de  Bruijn  [  h],  Bellman  and  Cooke  [2  ],  Hahn  [  5  ], 

Halanay  [  6  ]  and  Shimanov  [17,18])  The  author  has  been  influenced  by  all 
of  these  authors,  but  especially  by  Shimanov. 

Consider  the  equation 

0 

(3.1)  x(t)  =  /  [dT)(e)]x(t  +  0),  tiO 

-r 

and  its  formal  "adjoint" 

0  T 

(3.2)  y(s )  =  -  /  [dt|  (9)]y(s  -  e),  6  s  0 

-r 

T 

where  x,  y  are  n-vectors  and  the  symbol  B  always  denotes  tho  transpose 

of  a  matrix  B.  The  term  adjoint  is  Justified  by  the  following  observa- 

0 

tion.  If  we  let  Ix(t)  =  x(t)  -  /  fdT)(0)]x(t  +  0), 

0  T  ,  'r 

k*y(B)  =  y(8)  +  /  [in  (0)ly(s  -  9)  be  operators  defined  on  differentiable 

-r 

functions,  then 
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yT  (t)(Lx)(t)  +  [(L*y)(t)fx(t) 


|f[yT  (t)x(t)  -  /  /  yT  (t  +  |  -  0)[dri(0)]x(t  +  fc)dj 

~r  0 


The  expression  on  the  right  hand  side  of  this  equation  will  play  an  im¬ 
portant  role  in  the  following  discussion  and,  thus,  we  give  it  a  special 
designation.  For  any  9  €  C(  [-r,  0],  Rn),  t  €  C(  [0,  r],  Rn),  we  define 
the  symbol  (+,  9)  by  the  relation 


(30) 


T  0  0  m 

(♦>  9)  -  t  (0)q>(0)  -  /  J  *  (£-^)[dTi(0)l<p(i)dt. 

-r  0 


For  the  solution  x(q>)  of  (3.1)  with  initial  value  9  in 
C([-r,  0],  R11)  at  zero,  we  have  already  defined  the  operator  J(  t), 
t  2  0,  by  the  relation  xt(9)  =  ^ (t)<p  and  the  infinitesimal  generator 
{{  of  3 (t)  was 


(5.10 


<p(a)  =  -  r  s  0  i  0 

7  ♦(<>)-  /  [dT)(0)]f(0) 


If  t  is  in  C(  I0,  rj,  Rn),  then  system  (3*2)  has  a  solution  y(i|r) 

with  initial  function  +  at  zero  and  defined  for  all  s  £  r.  If  we 

let  yg(t)(e)^^y(y)(8  +  0),  0  *  e  *  r,  s  £  0,  then  the  operator 

»  s  0,  defined  by  yff)  *J/  *(b)1  has  all  of  the  same  pro- 

s 

perties  as  *7  (t).  The  infinitesimal  generator  of  is  de- 


dv>  -  11m  _  j  [7*(a)t  -  tl. 


8  -»0 


fined  by 
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For  later  purposes,  it  is  convenient  to  define  an  operator  $  *  on 
C(  [0,  r],  Rn)  in  such  a  way  that  A  direct  calculation 

shovs  that 


y(e),  o  s  e  s  r 

(3.5)  =  *1  -  o 


-♦(o)-/  [dij  (e)Jt(-e). 

-r 


It  is  also  easy  to  verify  that  v^*(s)  and  satisfy  the 

following  relationship: 


(3-6) 


B  £  Q. 


It  is  obvious  that  x(t),  y(t)  continuously  differentiable 
n-vector  functions  and  x^.,  y^  defined  by  x^(0)  =  x(t  +  0), 
yt(0)  =  y(t  -  0),  -  h  i  0  s  0,  implies 

dt^t ’  t 9  xt^  +  ^yt'  xt^ 


where  y^(0)  =  dy(t  -  0)/dt  =  y(t  -  0),  £^.(0)  =  dx(t  +  0)/d0  *  x(t  ♦  0), 
-  h  £  0  S  0. 


Lemma  HI.l, 


(♦,  <f9)  =  <P)  for  all  9  In  #«?>, 


♦  in 
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£22?:  (♦,  -  tT  (o)«Ko)  -  /°  /  tT  ({  -  e)[dti(e)  JKi)d| 

-r  0 

=  ♦*  (o)4»(o)  -  /°  [yT  (i  -  e)(dn(0)(pU)]J 

-r 

+  /°  /  dt"-ip  9}  [dT)(e)j<p(s)d£ 

-r  0  5 

0  0  8  T 

=  /  *T  (-e)(dti(e))<p(o)  +  /  /  dMf— ^'[dri(e)  >p(£)d£ 

-r  -r  0  5 

0  Q  T 

=  -  tT  (0)9(0)  -  /  /  [  -  d*.  If  -  ^]dn(e)<p(£)d£ 

-r  0  s 

=  (<1%,  <p) 
as  was  to  be  shown. 

Lemma  in.  2.  If  t  ^  0  is  any  given  constant  and  q>  e  C( [-r,  0],  RR), 

\|r  €  C([t,  t  +  r],  Rn),  then 

(f/*(t  -  r^(t)<p)  =  constant  for  0  5  t  $  t. 

Proof:  From  the  properties  of  and  v7,  it  follows  from  Lemma  IU.1 

that 

J^(f7*(t-T)\|r,  ^(t)<p)  =  -(^V/*(t-T)y,  +  (^(t-T)i|r,  ^J( t)<p) 


Lemma  HI.  3.  If  X  ^  u  then  for  any  nonnegative  integers  k, 
y  e  ({*  -  nl)k+1,  and  q>  e  iff{  -  Xl)l+1  implies 


(((T  -  ul)k_Py,  (<^  -  Xl)i_q<p)  =  0, 


0  S  p  *  k,  0  *  q  *  i. 
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Proof:  We  always  suppoee  X  ^  p.  Suppose  k  =  0.  We  wish  to  show  that 
that  for  any  nonnegative  integer  1  and  any  q,  0  £  q  £  I, 

(3-7)  (i,  (<^  -  Xl)l-q9)  =  0. 


We  prove  this  by  induction 
^9  =  Xcp  which  implies 


on  i. 


(*,  <^<p)  =  x(y,  9)  =  {({%,  q»)  =  u(*>  <p) 


and  thus  (y,  9)  =  0.  If  we  suppose  (5*7)  has  been  proven  for  0  *  1  £  T  -  1* 
then  for  y  -  pi),  ( -  Xl)r~S>  €  X)*(+1,  the  induction 

hypothesis  implies 


(3.8)  (♦,  (<#  -  Xl)r-q9)  =  0  for  0  *  q  *  r  -  1. 


Furthermore,  for  q  =  y  -  1, 

o  =  (♦,  U(  -  xi)q»)  =  -X(V,  9)  +  <p)  =  U  +  n)(t,  9) 

which  implies  (y,  9)  =  0.  But  this  is  (5*6)  for  q  =  f.  Consequently, 
relation  (5*7)  is  true. 

Nov  let  us  suppose  the  conditions  of  the  lemma  are  satisfied  for 
all  k  *  x  -  I  and  all  integers  /.  We  wish  to  show  this  is  true  for 
k  =  r*  Since  for  k  =  y  and  any  1,  y  €  ff  -  pl)k+1,  9  e  ^(^-Xl)l+1 

implies  -  ul)r"py  e^f(^*  -  pi)1*1,  (^-  XI)l-Q9  €  ftf(  ^  -Xl)q+1, 

for  0  -  p  -  y ,  0  S  q  5  1,  the  induction  hypothesis  implies  that 


(3.9)  ((#*  -  pl)r_Py,  (<7->>l)i"S))  =  0,  0«p<r-l,0MM. 
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It  remains  to  shew  these  relations  are  true  for  p  ~  y;  that  13, 

(3.10)  (t,  (^-XI)^)  =  0,  0  s  q  S  /. 

Let  us  first  show  that  (\|r,  (^-XI)*9)  =  0.  We  know  from  (3*9)  that 
0  =  ((4*  -  Wl)t,  (tf-  XI)1*)  =  (*,  (<^-ul)(  <^-Xl)V  =  (x^i)(t,  (d  -Xl)*q>) 

since  (^-Xl)*9  =  X(  -XI)  *9.  This  is  the  desired  relation.  Con- 
sequently,  (3*10)  is  true  for  q  =  0.  Suppose  (3*10)  has  been  shown  to 
he  true  for  all  q  £  t)  -  1.  From  (3.9)* 

0  =  (ttf  ill)*,  =  (♦,  (4* i  +  (x-mJDC^-xi)1-^)  = 

=  (♦,  (^-xi)i-(,l-1)9)  +  (X  -  M,  («^-x i)f-\) 

«(x-n)(t,  ( d(  -xi)i-\) 

by  the  induction  hypothesis.  Consequently,  (+,  (^-*J)*“\)  =  0  and 
(3.10)  is  valid  for  all  q.  This  completes  the  proof  of  the  lemma. 

Corollary  HI.l.  If  +  €  nfj.cC),  *«<(<?>.  then  (*,  9)  =  0 
if  X  /  >i. 

Proof:  Take  p  =  k,  q  =  f  in  the  above  lemma. 

Lemma  III  A.  If  t  €  then  for  all  € 


(*,  (  d  -«)»)  ■  Oi 


of 


that  Is,  ({(*  -  XI)  Is  orthogonal  to  the  range, 

^  -XI,  the  orthogonality  being  with  respect  to  the  symbol  (+,  9) . 
If  9  €  tfta-K),  then  for  all  +  € 

M*  -  «>*.  »)  -  0- 


Proof:  This  is  an  immediate  consequence  of  Lemma  III.l. 

Lemma  III. 5.  X  is  in  Pa(^)  if  a  only  if  X  is  in  Pcr(^* ) . 

Proof:  X  is  in  Po(^)  if  and  only  if  det  A(\)  =  0  where  A(x)  is 
defined  in  (2.6).  X  is  in  Pcr(^*)  if  and  only  if 

t(e)  =  eAeb 


and  b  is  a  nonzero  solution  of  the  equation. 


0 


[XI  -  /  [dT]1  (a)]eXs]b  =  A1  (x)b  =  0, 

~r 


or  if  and  only  if  det  AT  (x)  =  0.  But  det  AT  (x)  =  det  A(x)  and  the 
lemma  is  proved. 

Lemma  IU.l,  Corollary  III.l,  Lemma  III. k  and  Lemma  HI. 5  indicate 
that  the  operators  and  are  in  some  sense  adjoint  to  each  other; 

namely,  (*,  ({<»)  =  {&**,  <p),  ^  *  and  ^  have  the  same  point  spectra®, 

?ffY  (SC)  Is  orthogonal  to  if  X  /  p  and  C  -  XI)  is 

orthogonal  to  t\{  ^  -  XI).  If  we  could  show  that  wfxi<b  and 
w  m')  have  the  same  dimension  and  the  matrix  formed  by  (+j,  9^), 
where  the  i/ y  are  bases  for  Jjfyifl),  respectively,  is^thenlliere’ls 
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almost  a  complete  similarity  with  the  usual  concept  of  adjoint.  The 
next  results  are  a  precise  statement  of  these  remarks. 

Lemma  III. 6.  If  X  is  in  Pa(^),  then  the  dimensions  of  ^Yyk^L) 

and  are  equal,  -  Xl)k,  $fx(d)  =  <^-Xl)k, 

for  8ome  k  with  -  Xl)k  /  -  Xl)*"1, 

iff  (  4  -  Xl)k  ft  %f  (  -  Xl)k.  If  q>  «  then  a  necessary  and 

sufficient  condition  that  9  €  Xi)k  is  that 

(y,  <p)  =  0  for  all  y  e  {fl*  -  \l)k. 

Proof:  First,  we  introduce  some  notation.  With  the  matrix  A(\)  defined 
in  (2.6),  we  define  matrices  P^(x)  as 

pj+1(x)  =  ~j7^  '  a(J)(^)  =  ~y-  (AU))>  J  -  0,1,2,  ... 

and  the  matrix  of  dimension  (kn)x(kn)  as 


l£t  =  col  (d^,  ...,  a^),  0^  =  col  {0^,  ...,  0^),  where 

each  o^, 


is  an  n  x  p  matrix,  be  bases  for 


respectively. 


Suppose  9  €  TC{  c(-  \i)k.  Then  it  is  easy  to  show  that 


<p(e) 


k-l 


2  r 

J=o 


J+i 


J.9 

~ re 


inhere  A^r  *  0,  r  =  col  (r^,  . ..,  Tk)j  that  is,  the  vector 
belong  to  the  null  space  of  A^.  Consequently,  a  basis  9^ 
T{({(-  xi)k  can  be  written  as 


X  must 
for 


k=l 

Me)  =  2 

k  J=o 


Ct  ^  e-w 
T*,J+1  J|  e  ■ 


In  the  same  way,  one  showB  that  a  basis  for  -  Xl)^  ^an  be 

written  as 


v»> 


k-l 
=  Z 
j=0 


a  ^ 

%k-j  y. 


These  remarks  prove  that  the  dimensions  of  2^(  (X  -  xi)k  and  tftf-x  Dk 
are  the  same  for  every  value  of  k  and  the  first  part  of  Lemma  III. 6 
follows  immediately. 

Now,  suppose  that  q>  e  (  d-  XI)1  and  if  €  W-I,  I)k. 

Then  t(e)  =  Tk(e)h,  qp(e)  «  0f(e)a,  for  some  constant  vectors 
b,a,  and 


0  9 

(*,  ?)  =  t>T  TkT(0)*f(0)a  -  /  /  bT  YkT(i-0)d'1(e)*i(s)a  d{ 

■  >T  -  £  ‘t  (u[|>,,(,)  /* 


J=o 

k-l  i-1  T 


0  01+J+1 


=  b  [pkk“kl  '  iJ)  ^  ^_r  IT+3+1)-*  eX0dti(0))a^J+1]a. 
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Consequently, 


(3.11) 


<*>  »>  ■ bT  *  ■  v<  *  ■  V 


For  any  9  €  iff  {  -  XI)*,  v©  wish  to  derive  necessary  and  suffi¬ 

cient  conditions  that  9  €  td  -  XI )k;  that  is,  there  exists  a  9^  such 


that 


(3.12) 


( {{  -  XI)  V  =  <p- 


If  9  e  d^(  -  Xl)t  there  existe  an  a  such  that 


(3.13) 


<p(0)  =  *4(e)a  =  ai,  j+ia  Jr e  =  ^  rj+i  jr e  » 


rj+i =  ai,j+ia'  J  =  0,1>  i-1* 


If  (5.12)  is  to  have  a  solution,  then  necessarily 


(55  -  x)\(0)  =  9(e) 


or 


(3-14) 


n+ifre  ^iki).-eXe> 


where  the  are  n -vectors  satisfying  some  additional  conditions 

which  we  proceed  to  derive. 

If  we  define  9*m;(0)  to  be  equal  to  (d/d0  -  X)mq>#(0)  with 
9  given  in  (3.14),  then 
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<piB)(e) 


k-m-1 

jS> 


ek+J^  xe 


eX0  +  lr1  r  * _  -  - 

^  rj+l  (*+3^).'  a  ’ 


m  *  0, 1,2,  .  * .,  k  -  1. 


x«o 


If  is  to  satisfy  (3.12),  then  '  must  be  in (  0(  -  XI ) 
for  m  =  0,1, 2,  . ..,  k  -  l;  that  is, 


*<m)(0)  =  /  dT,(a)fW(e) 


-r 


or, 


+  r*  ~ 
'm+2 


0 

-/  e 

-r 


X0 


dT)(0)( 


k-m-1 

Z 

J=o 


„  0 

T5+J+1  J7 


I-1  gfc+J-m 

^  rj+i  TE*33).,)' 


m  =  0,1,2,  ...,  k  -  1. 

Using  the  definition  of  the  matrices  P  above,  these  last  equations  can 

J 

be  written  as 


1=1 

Plrm+1  +  P2rnH-2  +  •••  +  Purk  =  ^  Pk+j-nt+lrj+r 

Letting 

r*  =  coi  (rj*,  r{) 

d  =  col  (dj/  •••;  d^) 

1-1  |-1 

ds  =  j5o  Pjt+J-a+2rJ+1  =  £q  Pk+j-8+£?|,j+ia> 

we  hare  the  following  result:  a  necessary  and  sufficient  condition  that 
?  c  Tt  {.&s  XI)1  incites  <p  £  Xl)k  is  that  the  set  of  eqmtions 
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V* 


=  d 


have  a  solution;  that  is,  eT  d  =  0  for  all  elements  e  €  ).  But, 

for  any  e  e  (A^  ),  there  exists  a  b  such  that  e  =  0^b,  so  that 

T  T 

the  desired  property  holds  if  and  only  if  b  3^  d  =  0  for  all  p-vectors 
b,  or,  if  and  only  if, 


0  = 


A 
b  Pkd 


It  T  1-1 


b  jJj  Pk+j-8+2ai,J+l)a 

bT  (Bfi  j5)^  k-iPi+J+2ai,j+l)a 


for  all  b.  But,  from  (3*11);  this  latter  relation  is  true  if  and  only 
if  (*r,  q>)  =  0  for  all  *  €  7? {4*  -  Xl)k. 

This  completes  the  proof  of  the  lemma. 


lemma  III. 7.  If  f  is  a  basis  for  Aft*)  and  <J>  is  a  basis  for 
^),  then  (f,  $)  is  nonsingular  and  may  be  taken  as  the  identity. 
Furthermore,  if  (f,  $)  =  I  and  the  square  matrices  B*,  B  are  defined 
by  =  *B*,  #  &  =  4>B,  then  B*  =  B  T. 

Proof;  Suppose  «'*T  yii'C)  =  ~  kl)k,  )  =  J(  (rf  -  XI  )k 

where  k  is  the  least  integer  for  which  this  is  true  and  f,  <J>  are  bases 
for  WXU(*),  XU[),  respectively.  If  there  exists  a  constant  vector 

a  such  that  (f,  C)a  =  (f,  $a)  =  0,  then,  for  <p  =  <&a,  we  have 
(+,  cp)  =  0  for  all  Jr  €  t  C  (o/  *  XI)  .  Consequently,  from  lemma  III. 6, 

*  *  i’<  a.  XI )k.  Thus,  there  is  a  <p#  such  that  <p  =  -  Xl)^* 

Obviously,  and  “^(^-Xl)^  implies 
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9  =  ( /X  ~  M)k<P*  =  0.  Finally,  0  ■  (p  «  «a  implies  a  =.  0  since  % 
is  a  basis  for  W ^{(l)  and  the  matrix  (Y,  #)  is  nonsingular.  By  a 
change  of  basis  one  can  obviously  take  (Y,  #)  =  I.  From  lemma  III.  1 
and  the  definition  of  (Y,  9),  we  have 

(Y,  4*)  =(*,  *B)  =  (Y,  *)B  -  B 

-  (d\  •)  =  (fB*,  •)  =  (f,  *)  -  B*  T, 


which  cample  tea  the  proof  of  the  lemma. 

If  one  does  not  chooBe  (f,  ♦)  =  I,  then  the  matrices  B*  and 
B  are  related  by  B  =  (f,  $)"1B*T  (?,  #). 

Lemma  111.8.  If  X  is  In  Pa(/f)  and  $  Is  a  basis  for  *^(<^/),  then 
the  solution  x(qp)  of  (3*1)  with  initial  function  <p  in  at  zero, 

is  defined  for  all  values  of  t  in  (~®,  *)  and 

x(9)(t)  =  -  «  <  t  <  », 


p 

where  B  is  the  matrix  defined  by  *  tB.  If  f  is  a  basis  for 
yT\(c(*)’  then  the  8olution  y(t)  °£  (3«2)  with  initial 

function  \|r  in  ifl  defined  for  all  t  in  (-»,»)  and 

T 

y(+)(t)  =  Y(0)e-B  t,  -  “  <  t  <  ® 


where  B  is  the  matrix  defined  by 


c(*t 


TB 


The  proof  of  this  is  obvious. 

Let  i()>  «  ^(-4*)  b0  defined  as  before,  have  dimension 

p,  and  let  !  be  bases  of  respectively.  Pram 

Lemma  III. 7,  ^  can  choose  (f,  <&)  ~  I,  the  identity.  For  any 
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9  6  C(  [-r,  0],  Rn),  ve  define  a  rector  b  e  i?  and  a  function  (p  ;  c([ 
by  the  relation 

(3.15)  9  =  $b  +  9,  b  *  (f,  9). 

It  follows  that  (f,  f)  a  0  and,  also,  this  decomposition  is  unique. 

If  we  define  by  £^(0)  =  i(t  +  0),  -  r  ^  0  i  0,  then  the 

equation  (3-1)  can  be  written  as 

(3.16)  ^  t  10. 

For  any  solution  x^  *  2^(9),  9  111  C( [-r,  0],  Rn),  consider  the  change 

of  variables 

(3.17)  \  53  ®y(t)  +  2t,  t  S  0,  y(t)  =  (f,  xt). 

Since  (f,  z^)  =  0  for  all  t,  this  implies 

t(t)  =  |f  (t,  xt)  =  (f,  tfct)  *  B(r,  xt)  =  Ry(t) 

*t  =  *t  -  ^t(t)  »  -  ®3y(t)  -  ^(xt  -  *y(t))  =  C(zt, 

(^>  z^)  *  o. 

Consequently,  ve  have  the  following 

Theorem  m.l.  Suppose  the  operators  ,  //*  are  defined  by  (3.4), 

(3.5 )>  respectlTDly,  and  let  -Wx(4)»  -)YX(4*)>  Xe  Pc(^),  be  the 
naximal  subspaces  of  C( [-r,  0],  Rn),  C([0,  r],  Rn)  annihilated  by  powers 
of  (^  -  XI),  (<^*  -  XI),  respectively.  If  ♦  is  a  basis  for 
T  18  a  ba8lfl  f0r  (*>  ♦)  “ 


the  identity,  where 
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(♦>  <P)  1°  defined  in  (3.3)>  then  the  change  of  variables  (3.17)  applied 

to  (3.16)  yields  the  equivalent  system 


(3.18) 


y(t)  -  B  y(t) 

£(t)  «  (?,  zt)  -  0 


where  B  is  a  matrix  defined  by  =  $B,  ha 8  all  of  Its  eigenvalues 
equal  to  X,  and  the  point  spectrum  of  the  operator  restricted  to 
the  set  of  9  such  that  (f,  9)  »  0  does  no ,  contain  X. 

Now  let  us  consider  the  perturbed  equation 


i(t)  «  f(xt)  ♦  G(t,  xt), 

(3.19)  0 

f (9)  -  I  UTj(e)]9(e), 


where  G  is  Bone  function  defined  for  0  S  t  S  •,  9  €  C([-r,  0],  Rn). 
If  we  define  the  operator  l\  as  in  (3.4)  and  the  operator  by  the 

relation 


(3.20) 


-r  i  9  <  0, 

0-0, 


then  syBtem  (3*19)  is  equivalent  to 


(3.21) 


With  ♦,  f  defined  as  before,  (f,  ♦)  «  I,  the  transformation 


(3.17)  applied  to  (3.21)  yields 


-  32  - 


*(t)  -  (f,  ±t)  -  (?,  d\)  + 

-  ( 4  V  xt)  + 

‘  B(T,  xt)  +  (»,^txt) 

“  By(t)  +  (f,^xt), 

*t  “  *t  -  •S’(t)  -  4xt  +5;41It  -  *By(t)  -  *(f,^xt) 

*■  4(xt  “  *y(t))  +7t.xt ' 

=  ^2t  +^xt  -  *(^ 

But,  a  simple  calculation  shows  that 

(T,  4;xt)  +  *  1  xt> 

and  we  have  the  following  t'  eorem. 

Theorem  III, 2.  Suppose  the  operators  4.  A'  are  defined  by  (3.4),  (3«5)> 
respectively,  aid  let  »'(<?),  <<<**>,  x  €  Pa(^),  be  the  maximal 
subspaces  of  C([0,  r],  Rn)  annihilated  by  powers  of  (  -  XI),  {tf*  -  XI), 

respectively.  If  0  is  a  basis  for  T  ie  a  basis  for  <$/^(^*), 

(f,  9)  defined  in  (3.3)*  then  the  change  of  variables  (3. 17)  applied  to 
system  (3.20),  (3*21)  yields  the  equivalent  system 

y(t)  =  By(t)  +  fT  (0)G(t,  4>y(t)  +  zt), 

4t "  (lzt  +<4(®y(t)  f  zt)  "  «T(o)G(t,  *y(t)  +  zt)> 

(T,  zt)  =  0 


(3-22) 
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where  B  i a  the  matrix  defined  by  $  %  ■  #B,  has  all  of  its  eigenvalues 
equal  to  X,  and  the  point  spectrum  of  the  operator  C(  restricted  to  the 
set  of  9  such  that  (?,  9)  *  0  does  not  contain  X. 

Remark  UI,1.  By  a  repeated  application  of  the  above  process,  it  follows 
from  Corollary  m,l  that  one  make  a  further  decomposition  of  the  space 
C([-r,  0],  Rn)  which  will  yield  a  system  of  the  form  (3*22)  where  the 
real  parts  of  the  point  spectrum  of  the  operator  ^  restricted  to  the 
set  of  9  such  that  (l,  9)  -  0  are  less  than  any  preas signed  value, 
p.  The  matrix  B  will  then  have  its  eigenvalues  equal  to  the  elements 
of  P tj(0)  which  have  a  real  part  &  p. 


IV.  Perturbation  of  linear  systems.  In  this  section,  v*  indicate  some 


applications  of  the  results  of  the  previous  section  to  the  system 

0 

(b.l)  ±(t)  =  /  [dT,(®)Jr(t  +  9)  +g(t) 

-r 

where  x  is  an  n -vector,  t \(q)  is  an  n  x  n  matrix  whose  elements  are 
of  bounded  variation  on  [-r,  0]  and  g(t)  is  continuous  on  (-»,»). 
Together  with  system  (4.1),  we  consider  the  adjoint  equation 

(4.2)  *(t)  -  -  J°  [dT,T(0)]y(t  -  e). 

-r 

Lena  IV.  1.  If  tiO,  .^T(o)  =  I,  is  a  strongly  continuous  semi¬ 

group  of  operators  of  a  Banach  space  into  Itself  and  if  for  some  r  >  0, 

the  spectral  radius  p  =  is  <  1  f  0,  then 

for  any  function  h  mapping  (-»,  «)  into  such  that  h(t)  is  almost 

periodic,  the  function 

(4.3)  Z*  =  /  (7(t  -T)h(T)dT 

-c* 

is  almost  periodic  in  t  with  the  same  frequencies  as  h(t), 

||z|||  £  KR/p,  -  20r  *  log  p,  K  constant,  R  =  supt)|h(t)  ||,  and  is  a  uni¬ 
formly  asymptotically  stable  solution  of 

(4.4)  *  C(  +  h(t) 

where  is  the  infinitesimal  generator  of  f/(t). 

Proof:  If  -2gr  *  log  p,  and  r  *  P,  then  Theorem  II.  1  implies 
||£/(t  -  t  )qp  J1  *  KB-P(t-T)(|<p(|  for  all  til,  and  some  constant  K.  There¬ 
fore,  the  integral  in  (4.5)  exists  and  ||z*||  %  KR/p.  Furthermore, 
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^  (0)  =  I  and  d<7(t)/dt  =  $3  (t)  Implies  z|  satisfies  (4.4).  The 

uniform  asymptotic  stability  of  z*  follows  from  Theorem  II.  1.  To  show 

that  z*  is  almost  periodic  with  the  same  frequencies  as  h(t),  it  is 

sufficient  to  show  that  for  every  sequence  of  real  numbers  (t  )  for 

m 

which 

h(t  +  t^)  -  h(t)  -*  0  as  m  ->  uniformly  for  t  in  {-*>,  «), 
we  have 


z*_^  -  z^  ->0  as  m  -*«,  uniformly  for  t  in  »). 

m 

This  is  easily  verified  using  Theorem  II.  1  and  the  lemma  is  proved. 

The  proof  of  the  next  two  lemmas  are  standard  and  may  be  found  in 

[7,12]. 

Lemma  IV.2.  If  C  is  a  constant  n  x  n  matrix  whose  eigenvalues  have 

real  parts  £  2p  >  0,  and  h(t)  is  an  almost  periodic  n-vector,  then 
there  is  a  unique  almost  periodic  solution  of  the  equation 

w  »  Bw  +  h(t) 

which  is  bounded  by  KR/p  for  some  constant  K  and  R  =  supt|h(t)|. 

Lemma  IV.  If  B  is  a  constant  n  x  n  matrix  whose  eigenvalues  may  have 

zero  real  parts,  h(t)  is  an  almost  periodic  n-vector  with  a  finite 
Fourier  series,  then  a  necessary  and  sufficient  condition  that  the  syBtem 


f  =  By  +  h(t) 
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have  an  almost  periodic  solution  is  that 

1  T 

11m  s  /  *.(t)h(t)dt  =  0,  J  =  1,2,  ...,  r, 

T  -+«•  0  J 

where  ...,  tr)  1b  a  basis  for  the  almost  periodic  solutions  of 

the  adjoint  equation 

-  yB 

One  could  also  state  a  result  similar  to  IV.  3  for  the  case 

when  h(t)  does  not  have  a  finite  Fourier  series,  hut  it  is  necessary  to 
have  an  hypothesis  which  guarantees  that  the  integrals  of  certain  almost 
periodic  functions  are  almost  periodic.  We  do  not  discuss  this  question 
here. 


Lemmas  IV.  1  -  IV.  3  together  with  Theorem  III.l  and  Remark  III.l 
yield  in  a  very  natural  way  extensions  to  functional-differential  equations 
of  the  standard  results  concerning  the  existence  of  periodic  and  almost 
periodic  solutions  of  ordinary  differential  equations  which  are  perturba¬ 
tions  of  &  linear  system  with  constant  coefficients.  We  do  not  state  all 
of  these  results  for  functional-differential  equations  but  merely  give  some 
indications  of  the  manner  in  which  they  are  obtained.  After  inspection  of 
the  proof 8,  one  will  see  that  the  arguments  are  essentially  the  same  as 
for  ordinary  differential  equations. 

Theorem  4.1.  If  g(t)  is  27r-periodic  and  (f^,  ...,  represents  a 

basis  for  the  2rr-periodic  solutions  of  the  adjoint  equation  (4.2),  then 
a  necessary  and  sufficient  condition  that  system  (4.1)  has  a  2jr-periodic 


solution  is  that 
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/  K  (t)g(t)dt  «  0. 

0 

If  there  are  no  2ir-per  iodic  solutions  of  the  adjoint  equation,  then 
system  (4.1)  has  a  unique  2r-periodic  solution  bounded  by  KR/p  for 
some  constant  K,  R  »  8up^|g(t)|  and  2p  being  defined  as  the  largest 
number  such  that  |x  -  in|  2p,  n  =  0,  +1,  +2,  for  all  roots 

X  of  the  characteristic  equation  (2.6). 

Proofs  Let  (%*  be  the  operators  defined  by  (3.4),  (3*5).  l£t 

$1,  Y^  be  bases  far  the  linear  extensions  of  the  manifolds 

Turfjt/t'),  respectively,  for  all  X  in  Pcr(^)  with  X  ft  0  (mod  i),  ReX  *  0, 
and  let  #2,  Y^  be  bases  for  the  linear  extensions  of  the  manifolds 

yfyid*)*  respectively,  for  all  X  in  Po(^)  with  XsO  (mod  i). 
Furthermore,  we  can  choose  (Y^,  *1)  =  I,  (*2>  ®2)  “I,  "tere  (t,  9) 
is  defined  in  (3*3)  (see  Lemma  UI.7).  Consequently,  by  Theorem  III. 2 
and  Remark  III.  1  the  transformation  of  variables 

xt  =  ®iw(t)  +  +  zt» 

(U.5) 

(*!_»  \)  =  w(t),  (T2,  xt)  =  y(t), 

applied  to  (4.1)  yields  the  equivalent  syBtem 

fi-(t)  =  B^t)  +  Tir(0)g(t) 

(4.6)  fr(t)  =  ^y(t)  +  Yg  (o)g(t) 

£(t)  =  ((zt  -  VlT  (O)g(t)  -  *2T2T  (O)g(t), 

((f^), zt)  =  0, 
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where  ,^(0)  =  0  for  -  r  M  <  0,  *  g(t)  for  $  -  0,  1  ,  «  are 

the  matrices  defined  by  4  =  ^B,  4$2  =  respectively,  all 

eigenvalues  of  are  ^  0  (mod  i),  all  eigenvalues  of  are 
2  0  (mod  i)  and  all  elements  in  the  point  spectrum  of  the  restriction 
of  to  the  set  of  all  9  such  that  ((f^  9)  -  0  have  negative 

real  parts.  System  (4.1)  will  have  a  2ir-periodlc  solution  if  and  only 

if  there  is  a  2r -periodic  solution  w(t),  y(t),  z t  of  (4.6).  After 
observing  that  the  basis  . ..,  the  2?r-per  iodic  solutions  of 

the  adjoint  equation  (4.2)  are  of  the  form  f  ▼j(t)  where  the  v^(t) 

T 

are  2r-periodie  solutions  of  the  equation  f  =  -  B  v,  one  can  apply 
Lemmas  IV.1-IV.3  to  complete  the  proof  of  Theorem  IV.  1. 

In  the  same  manner,  one  proves 

Theorem  IV.  2.  If  g(t)  is  almost  periodic  in  t  with  a  finite  Fourier 

series  and  (9^,  ...,  f^)  represents  a  basis  for  the  almost  periodic 

solutions  of  the  adjoint  equation  (4.2),  then  a  necessary  and  sufficient 
condition  that  system  (4.1)  has  an  almost  periodic  solution  is  that 

lim  £  (t)g(t)dt  =  0,  J  *  l,2,...,k. 

T  T  0 

If  there  is  no  almost  periodic  solution  of  the  adjoint  equation,  then  there 
is  a  unique  almost  periodic  solution  of  (4.1)  and  it  is  bounded  by  KR/p 
for  some  constant  K,  R  *  sup^|g(t)  j,  and  2$  defined  as  the  largest  number 
such  that  |Re  x|  *  2p  for  all  roots  of  X  of  the  characteristic  equation 
(2.6). 

Ho  consider  the  nonlinear  eqmtion 

0 

(4.7)  *(t)  =  /  [dr)(e)]x(t  +  0)  +  G(t,  e) 

-r 
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where  i\($)  is  the  same  matrix  as  before,  €  is  &  parameter,  q(t,  <p,  <») 

Is  continuous  in  t,  q>,  €  for  -  »  <  t  <  »,  q>  e  C( [-r,  0],  Rn), 

Ml  -  H,  H  >  0,  0  ^  ^  cq,  and  is  Lipschitzian  in  9.  Furthermore, 

there  exists  a  function  q(€,  p),  continuous  for  Q  £  €  -  €q,  0  *  p  s  h, 
such  that  tj(0,  0)  =  0  and 

|G(t,  <p,  e)  -  G(t,  *r,  e)|  i  tj(c,  p) |j<p  -  til,  q(t,  0,  0)  *  0, 

for  all  q>,  t  c  C(  [-r,  0],  Rn),  M,  ||*r||  i  p,  0  £  €  s  €q,  -  <  t  < 

0 

If  the  characteristic  equation  (2.6)  of  k(t)  =  /  [dq($)]x(t  +  9) 

-r 

has  no  roots  =  0  (mod  i)  and  G  is  2rr-per  iodic  in  t,  then  one  can 
prove  that  there  is  a  unique  2rr-period5c  solution  of  (h.7)  in  a  neighbor¬ 
hood  of  the  origin  and  this  solution  approaches  zero  as  e  -►  0.  If  the 
roots  of  this  characteristic  equation  have  nonzero  real  parts  and 
G(t,  <p,  € )  is  almost  periodic  in  t  uniformly  with  respect  to  <p  for 
each  fixed  €,  then  there  is  a  unique  almost  periodic  solution  of  (k.T) 
in  a  neighborhood  of  the  origin  and  this  solution  approaches  zero  as 
€-4  0.  The  proofs  may  be  supplied  exactly  as  in  [  7;  Chapters  5,  12]  by 
making  use  of  the  change  of  variables  as  in  Theorem  17.1.  These  results 
have  been  obtained  previously  In  a  slightly  less  general  foim  (see,  for 
example,  Krasovskii  [10],  ghimanov  [16]). 

What  happens  in  the  case  when  some  of  the  roots  of  the  characteristic 
equation  (2.6)  are  s  0  (mod  i)  and  G  is  2r-periodic  in  t?  In  this 
case,  one  can  extend  the  method  of  casting  out  secular  terms  to  derive 
the  determining  equations  (bifurcation  equations)  associated  with  (h.J). 
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This  procedure  has  already  been  indicated  by  Shimancv  [  ],  but  ve 

repeat  it  here  for  completeness.  For  simplicity  in  notation,  ve  restrict 
ourselves  to  the  following  case:  if 

A(X)  -  XI  -  f  [dn(0)]eX0 
-r 

and  X  is  the  root  of  multiplicity  p  of  det  A(x)  =  0  which  i 3 
se  0  (mod  i),  then  there  are  n  linearly  independent  solutions  of 
A(x)*>  *  0.  Furthermore,  for  simplicity  only,  suppose  G  =  €H. 

5y  the  change  of  variables  (4.5),  system  (4.J)  can  be  written  as 

w(t)  =  BjV(t)  +  €  (0)0^ 

(4.8)  J’(t)  =  B^y(t)  +  e  f^(o)Gv 

\  +  e  -  e  *i?1T(0)G1  -  e 

((?!,  *2)\)  =  0, 

idle  re 

°i =  "W*  y^>  V =  +  +  zt)> 

{'k-9)  At(e) “  (  °  -  -rs8<°- 

I  ^  9  0  =  0, 

no  eigenvalue  of  is  m  0  (mod  i),  e  is  2>r -per  iodic  and  all 

elements  in  the  point  spectrum  of  the  restriction  of  to  the  set  of 
all  ^  such  that  ((f^,  F^),  qp)  =  0  have  negative  real  parte. 
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If  we  let  y(t)  *  e  v(t),  system  (4.8)  can  be  written  as 


Kt) » ^(t)  + 

(4.10)  +(t)  =  €e  **  1*(0)Gr 

4(t)  =  +€>lt  -  C«1f1TtO)G1  -  €4»2T2T(0)G1, 

((Tr  T2),  zt)  =  0, 

-But 

where  G^  =  G^(t,  w(t),  e  v(t),  z^).  One  can  now  use  Lemma  IV. 1  to 
repeat  all  of  the  arguments  in  [  7  ,  Chapter  6]  to  obtain  a  periodic  func 
tion  (v*(t),  v*(t),  z*)  with  a  vector  a,  the  mean  value  of  v*,  arbi 
trary,  and  a  set  of  equations  (the  determining  equations)  involving  the 
arbitrary  vector  a  which  have  the  property  that  (w*(t),  v*(t),  z*) 
will  be  a  periodic  solution  of  (4.10)  is  and  only  if  the  vector  a 
satisfies  the  determining  equations.  We  do  not  discuss  this  question 
further  since  all  details  are  easily  supplied.  This  procedure  extends 
the  method  of  Ce  sari  -Hale  -Gambill  to  functional -differential  equations 
of  the  above  type. 

In  some  problems,  It  is  more  convenient  to  introduce  polar  coor¬ 
dinates  for  the  vector  y  in  (4.8).  In  the  following,  we  may  allow  G 
in  (4.2)  to  be  almost  periodic  in  t  and  need  only  assume  that  the 
matrix  has  the  form 


(4.11) 


=  dlag  (A^,  • A^) 


k 


where  each  <jj  is  positive.  If  /  «  (y^,  ....  y^),  the  transforaa- 
tion 

y2j-l  =  PJ  eln  &J 

(4. 12)  jgj  *  Pj  cob  0jf  J  *  1>2j  ...t  if 

Pj_  *  (Pj_>  •••>  Pj^)>  ®  =  (Sjy  ®jj.) 

applied  to  (4.8)  yields  a  set  of  equations  of  the  fora 

6(t)  -  d  +  €  0(t,  0(t),  p(t),  w(t),  zt,e) 
p(t)  »  c  R(t,  0(t),  p(t),  w(t),  zt,  € ) 

(4.15)  *(t)  -  BjV(t)  +  €  V(t,  e(t),  p(t),  v(t),  zt,  e) 

4t  *  ^zt  +  e  z(%*  V  e) 

where  zt  is  required  to  satisfy  ((f^,  f2),  z^)  =  0,  and  d  =  (cr^ 
Equations  (4.15)  are  of  the  same  type  that  hare  been  considered  in  [  3,7  ] 
for  ordinary  differential  equations.  It  is  a  simple  matter  to  extend  the 
method  of  averaging  in  [  3,7  ]  to  systems  of  the  form  (4.15)  and  to  prove 
the  existence  of  Integral  manifolds  which  are  generated  by  equilibrium 
points  of  the  averaged  equations.  We  do  not  go  in  detail  on  these  general 
questions,  since  these  extensions  will  be  clear  to  the  reader  who  is 
familiar  with  the  results  for  ordinary  differential  equations. 
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